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Abstract:  Eligher-order,  large-domain  electromagnetic  simulations  are  presented  based  on  the  finite 
element  method  (FEM),  method  of  moments  (MoM),  and  physical  optics  (PO).  The  simulations  combine 
higher  order  geometrical  modeling  and  higher  order  field/current  modeling,  which  is  referred  to  as 
double-higher-order  modeling.  The  examples  demonstrate  that  accurate  and  efficient  FEM,  MoM,  and  PO 
simulations  require  both  higher-order  geometrical  flexibility  for  curvature  modeling  and  higher-order 
basis  functions  for  field/current  modeling  in  the  same  method.  It  is  optimal  to  have  the  geometrical  orders 
and  current/field  approximation  orders  of  the  elements  entirely  independent  from  each  other,  so  that  the 
two  sets  of  parameters  can  be  combined  independently  in  a  double-higher-order  model. 

Key  words:  electromagnetic  analysis,  finite-element  methods,  moment  methods,  physical  optics,  higher- 
order  modeling. 


1.  Introduction 

This  paper  discusses  higher  order  modeling  of  geometry  and  higher  order  modeling  of  fields  and  currents 
in  the  framework  of  the  finite  element  method  (FEM),  method  of  moments  (MoM),  and  physical  optics 
(PO).  Traditional  FEM,  MoM,  PO,  and  hybrid  simulation  tools  are  low-order  or  small-domain 
(subdomain)  computational  techniques.  In  the  small-domain  approach,  structures  are  modeled  by  volume 
and  surface  geometrical  elements  (finite  and  boundary  elements)  that  are  electrically  very  small 
(typically,  on  the  order  of  A./1 0  in  each  dimension,  X  being  the  wavelength  in  the  medium)  and  the  fields 
and  currents  over  the  elements  are  approximated  by  low-order  (zeroth-order  and  first-order)  basis 
functions.  This  results  in  a  very  large  number  of  unknowns  (unknown  field-  or  current-distribution 
coefficients)  needed  to  obtain  results  of  satisfactory  accuracy,  with  all  the  associated  problems  and 
enormous  requirements  in  computational  resources.  In  addition,  commonly  used  finite  and  boundary 
elements  are  in  the  form  of  cells  (bricks,  tetrahedra,  and  triangular  prisms)  with  planar  sides  and  flat 
(triangular  and  quadrilateral)  patches,  and  thus  they  do  not  provide  enough  flexibility  and  efficiency  in 
modeling  of  structures  with  pronounced  curvature.  According  to  the  higher-order  or  large-domain 
computational  approach,  on  the  other  side,  a  structure  is  approximated  by  a  number  of  as  large  as  possible 
geometrical  elements,  and  the  approximation  of  field/current  components  within  individual  elements  is  in 
the  form  of  a  single  (three-  or  two-fold)  functional  series  of  sufficiently  high  order  (e.g.,  [1-5]).  This 
approach  can  greatly  reduce  the  number  of  unknowns  for  a  given  problem  and  enhance  further  the 
accuracy  and  efficiency  of  the  analysis  in  practical  applications. 

We  present  our  higher-order,  large-domain  FEM  [4],  MoM  [5],  and  hybrid  MoM-PO  techniques,  which 
are  referred  to  as  double-higher-order  methods  because  they  combine  higher  order  geometrical  modeling 
and  higher  order  field/current  modeling.  They  enable  using  large  curved  FEM  hexahedral  volume 
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elements  and  MoM  and  PO  quadrilateral  surface  elements  that  are  on  the  order  of  X  (e.g.,  1 .5A.  -  2  k)  in 
each  dimension  as  building  blocks  for  modeling  of  the  electromagnetic  structure  (i.e.,  the  finite/boundary 
elements  can  be  by  three/two  orders  of  magnitude  larger  in  volume/area  than  traditional  low-order 
elements).  Element  orders  in  the  model,  however,  can  also  be  low,  so  that  the  lower-order  modeling 
approach  is  actually  included  in  the  higher-order  modeling.  The  examples  demonstrate  that  for  FEM, 
MoM,  and  PO  modeling  of  general  structures  that  may  possess  arbitrary  curvature  it  is  essential  to  have 
both  higher-order  geometrical  flexibility  for  curvature  modeling  and  higher-order  basis  functions  for 
field/current  modeling  in  the  same  method.  In  addition,  it  is  optimal  to  have  the  geometrical  orders  and 
field/current  approximation  orders  of  the  elements  entirely  independent  from  each  other,  so  that  the  two 
sets  of  parameters  of  the  double -higher-order  model  can  be  combined  independently  for  the  best  overall 
performance  of  the  method. 


2.  Higher  Order  Volume  and  Surface  Geometrical  Modeling 


As  basic  building  blocks  for  volume  geometrical  modeling  in  FEM,  we  use  generalized  curved  parametric 
hexahedra  [4]  of  higher  (theoretically  arbitrary)  geometrical  orders.  A  generalized  hexahedron  (Fig.  1)  is 
analytically  described  as 


m  ‘^u  ^ v  “-w 

K«>  v,  w)  =  £  rXKr  («,  A  w)  =  2  X  2 


m  n  l 

U  V  w 


m=0  n= 0  /=0 


-\<u,v,w<\,  M  =  (Ku+l)(Kv+l)(Kw  +  \),  (1) 


where  rhr2,  ...,  rM  are  the  position  vectors  of  the  interpolation  nodes,  and  Km  K„  and  Kw  (Ku,  Kv,  Kw  >  1) 
are  geometrical  orders  of  the  element  along  «-,  v-,  and  w-  parametric  coordinates,  respectively.  The 
functions  Lf(u,v,w )  are  defined  as 

LKr  (u,  v,  w)  =  L £  {u)LK:  (v)LKr  (vv) ,  i  =  m  +  n{Ku  + 1)  +  l(Ku  +  \)(Kv  + 1)  + 1 ,  (2) 

where  LKm‘ ,  Lknr ,  and  Z,fw  represent  Lagrange  interpolating  polynomials. 


Fig.  1.  A  generalized  curved  parametric 
hexahedron  of  arbitrary  geometrical  orders  for 
FEM  volumetric  modeling. 


Fig.  2.  A  generalized  curved  parametric 
quadrilateral  of  arbitrary  geometrical  orders  for 
modeling  of  MoM  and  PO  surfaces. 
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Surface  geometrical  modeling  in  MoM  and  PO  is  carried  out  using  generalized  curved  parametric 
quadrilaterals  of  higher  (theoretically  arbitrary)  geometrical  orders  Ku  and  Kv  [5],  A  generalized 
quadrilateral  (Fig.  2)  is  analytically  described  by  the  2-D  version  of  Eqs.(l)  and  (2). 


3.  Higher  Order  Field  and  Current  Modeling 


The  electric  fields  inside  the  FEM  hexahedra  are  represented  as  [4] 

JV„-1  Nv  Nw  Nu  Nv- 1  Nw  Nu  N '  JV^-1 

^ Hijk  f  'lijk  ^vijk  f  vijk  ^  wijk  f  wijk 


i=  0  /=0  A=0 


(=0  7=0  A=0 


i=0  7=0  A=0 


(3) 


where  N,„  N„  and  Nw  are  the  adopted  degrees  of  the  field  approximation,  which  are  entirely  independent 
from  the  element  geometrical  orders,  Km  Kv,  and  Kw,  and  /  are  curl-conforming  hierarchical  vector  basis 
functions  defined  as 
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J  J  J  ou  dv  dw 

In  MoM,  electric  and  magnetic  surface  current  density  vectors  over  every  generalized  quadrilateral  in  the 
model  are  represented  as  [5] 


/L  ^  ®wjr  fuij  ^  ^  ^vij  fvij  ’  ^  ^  PUIJ  ^  ^  PVIJ  /v. 


(5) 


i=0  y=0  i=0  7=0  i=0  7=0  i=0  7=0 

where  ;V„  and  /V,.  are  the  adopted  degrees  of  the  current  approximation  and  /  are  divergence-conforming 
hierarchical  vector  basis  functions  defined  as 


f  = 

J  uij 


P(u)vJ 

J 


u'P(v ) 

f  = - - a 

J  vij  j  v 


J  =1  a  Xfl, 


(6) 


In  PO,  basis  functions  for  approximating  J%  are  divergence-conforming  two-dimensional  interpolatory 
Chebyshev-type  polynomials  given  by 


=  C,0)Tjv)  _Cv,TM,(fl)T 

j  mod  u  5  vij  j  mod 

J  u-ui  v  —  V  •  J  u-ui  v  —  V  • 


(7) 


where  T  mod  and  T  are  the  modified  and  regular  Chebyshev  polynomials,  respectively,  and  C,,  are  the 
normalization  factors.  The  zeros  of  modified  polynomials,  w(mod,  are  obtained  by  scaling  the  zeros  of 
regular  Chebyshev  polynomials  by  a  factor  of  cos{  7i/[2(/V„+l )]  },  and  similarly  for  v""“' . 

In  FEM  and  MoM,  the  same  functions  /  are  used  for  testing,  which  gives  rise  to  the  Galerkin  method.  In 
PO,  a  modified  point-matching  technique  is  applied  at  the  interpolation  points  of  the  basis  functions, 
which  makes  the  extremely  large  PO-PO  projection  matrix  in  the  general  MoM-PO  matrix  system  be  an 
identity  matrix,  and  thus  tremendously  reduces  the  computational  costs  associated  with  electromagnetic 
interactions  in  the  PO  region.  In  addition,  the  implemented  point-matching  testing  in  the  PO  region  has 
proved  to  be  accurate  enough  in  all  applications  and  is  faster  than  the  Galerkin  procedure. 
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4.  Numerical  Results  and  Discussion 

As  an  example  of  higher  order  FEM  modeling  of  curved  structures,  Fig.  3  shows  the  average  percentage 
error  in  calculating  ko  for  the  first  eleven  modes  of  an  air-filled  metallic  spherical  cavity,  versus  the 
number  of  unknowns  for  the  solution  (the  error  is  given  with  respect  to  the  exact  result).  The  numerical 
results  are  shown  for  a  higher  order  FEM  solution  with  the  sphere  modeled  by  a  single  generalized 
hexahedron  of  the  2nd  geometrical  order  (Ku  =  Kv  =  Kw  =  2)  and  two  higher  order  single-element  solutions 
using  a  generalized  hexahedron  of  the  4th  geometrical  order  (K„  =  Kv  =  Kw  =  4).  In  the  first  4th-order 
geometrical  model,  the  control  (interpolating)  points  that  do  not  belong  to  the  cavity  surface  (hexahedron 
faces)  are  arranged  to  define  an  inscribed  half-radius  sphere,  while  in  the  second  4th-order  geometrical 
model  the  inner  control  points  define  an  inscribed  cube  with  the  spatial  diagonal  equal  in  length  to  the 
cavity  radius.  In  all  three  models,  the  field-approximation  orders  are  varied  from  Nu  =  Nv  =  Nw  =  3  to  Nu 
=  Nv  =  Nw  =  7  (^-refinement).  Note  that  these  are  literally  entire-domain  FEM  models  (an  entire 
computational  domain  is  represented  by  a  single  finite  element).  We  observe  from  the  figure  excellent 
convergence  properties  of  all  higher  order  solutions,  namely,  that  the  computation  error  decreases  rapidly 
and  monotonically  when  using  ^-refinement.  We  also  observe  a  significant  additional  improvement  in 
accuracy  as  a  result  of  using  geometrical  modeling  of  the  4th  order  instead  of  the  2nd  order  geometrical 
modeling.  In  other  words,  it  is  impossible  to  /i-rcfinc  the  higher  order  model  with  the  2nd  geometrical 
order  below  about  1%  error  in  calculating  ko  due  to  the  inherent  geometrical  error  of  the  model,  whereas 
the  ^-refinement  in  the  model  with  the  4th  geometrical  order  brings  the  analysis  error  quickly  down  to  a 
fraction  of  a  percent.  Finally,  we  note  that  the  4th  order  geometrical  model  with  the  inner  points  on  a  cube, 
having  less  distorted  coordinate  lines  throughout  its  volume  and  thus  being  able  to  represent  the  fields 
more  accurately,  yields  by  an  order  of  magnitude  smaller  average  error  in  the  low-error  region  than  the 
model  with  the  inner  points  placed  on  a  sphere,  which  is  an  indication  of  the  importance  of  careful 
meshing  considerations  when  elements  of  higher  geometrical  orders  are  used. 

As  an  example  of  higher  order  MoM  modeling  of  curved  structures,  consider  a  spherical  dielectric 
scatterer  1  m  in  radius  in  the  frequency  range  10  -  600  MFIz  [5],  The  relative  permittivity  of  the  dielectric 
is  £r  =  2.25  (polyethylene).  Shown  in  Fig.  4  is  the  RCS  of  the  sphere  calculated  using  two  higher  order 
models  with  the  sphere  surface  being  approximated  by  means  of  384  generalized  quadrilaterals  of  the  first 
geometrical  order  (Ku  =  Kv  =  1)  and  6  generalized  quadrilaterals  of  the  fourth  geometrical  order 
(KU  =  KV  =  4),  respectively,  along  with  the  analytical  solution  in  the  form  of  Mie’s  series.  The  adopted 
electric  and  magnetic  current  approximation  orders  in  the  first  model  are  Nu  =  Nv  =  2  in  all  elements,  and 
the  resulting  total  number  of  unknowns  is  6144.  In  the  second  model,  Nu  =  Nv  =  5  in  all  elements,  which 
corresponds  to  a  total  of  2400  unknowns.  We  observe  a  very  good  agreement  of  the  numerical  results 
obtained  by  the  higher  order  model  with  first-order  geometrical  modeling  and  the  analytical  results,  with 
the  RCS  numerical  predictions  being  slightly  shifted  toward  higher  frequencies.  We  also  note  that,  even 
though  this  is  an  almost  small-domain  application  of  our  large-domain  MoM,  where  a  relatively  large 
number  (384)  of  elements  (with  relatively  low  current  approximation  orders)  is  needed  for  the  sphere 
surface  to  be  geometrically  accurately  represented  by  parametric  surfaces  of  the  first  geometrical  order, 
the  largest  quadrilateral  elements  in  the  model  are  0.391free-space  or  0.582diei  on  a  side  (2diei  =  2free-Space/^£r)  at 
the  highest  frequency  considered,  which  is  still  considerably  above  the  usual  small-domain  limit  of  0.1X. 
On  the  other  hand,  the  agreement  of  the  model  with  fourth-order  geometrical  modeling  with  the  exact 
solution  is  excellent  in  the  entire  frequency  range  considered.  Note  that  the  curved  quadrilateral  elements 
in  this  model  are  about  1 .572rrcc_space  or  2.352diei  across  at  the  highest  frequency.  Note  also  that  with  the  use 
of  the  model  with  fourth-order  geometry  modeling  and  fifth-order  current  modeling  the  number  of 
unknowns  is  reduced  about  2.5  times  and  the  accuracy  is  increased  as  compared  to  the  model  with  first- 
order  geometry  modeling  and  second-order  current  modeling. 
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Fig.  3.  Comparison  of  three  higher  order  single-  Fig.  4.  Normalized  radar  cross-section  [RCS/(a27t)] 
element  (entire-domain)  FEM  solutions  -  average  of  a  dielectric  (£r  =  2.25)  sphere,  for  two  higher- 
errors  for  the  first  eleven  modes,  k0,  of  a  spherical  order,  large-domain  MoM  models,  along  with  the 
cavity  against  the  number  of  unknowns.  exact  solution  (Mie’s  series). 


As  an  example  of  higher  order  PO  modeling  of  curved  structures,  consider  an  antenna  system  consisting 
of  an  array  of  nine  A/2  dipoles  in  front  of  a  circular  metallic  cylinder  (Fig.  5).  The  radius  of  the  cylinder  is 
7.5A.  and  its  height  1 2A,  The  dipoles  are  parallel  to  the  cylinder  axis  and  are  situated  1.25A,  from  the 
cylinder  surface  with  an  angular  separation  of  3.75  degrees  between  the  adjacent  dipoles  in  the  array.  All 
of  the  dipoles  are  center-fed  by  point  generators  of  equal  amplitudes  and  phases.  The  lateral  surface 
(barrel)  of  the  cylinder  is  approximated  using  144  quadrilaterals  of  the  second  geometrical  orders 
( Ku  =  Kv  =  2)  and  each  of  the  cylinder  bases  (caps)  is  represented  by  a  total  of  24  second-order  and  36 
first-order  ( K„  =  Kv=  1 )  geometrical  elements,  with  the  seventh-order  current  approximation  in  both 
parametric  coordinates  ( N„  =  Nv  =  7)  for  all  of  the  quadrilaterals  in  the  model.  All  quadrilaterals  are 
approximately  2A,  on  a  side,  which  is  20  times  the  usual  low-order  limit  of  A/10.  Each  dipole  is  modeled 
using  two  straight  wire  segments  with  fourth-order  current  approximations  ( Nu  =  4).  In  the  hybrid  MoM- 
PO  analysis,  the  dipole  array  is  in  the  MoM  region  and  the  cylinder  is  in  the  PO  region.  Fig.  6  shows  the 
computed  radiation  pattern  of  the  array  in  the  horizontal  plane.  The  results  obtained  by  the  higher  order 
MoM-PO  technique  are  compared  with  the  results  of  the  pure  MoM  analysis.  We  observe  an  excellent 
agreement  of  the  two  sets  of  results  in  the  front  region  (for  angles  up  to  90  degrees).  In  the  back  region 
(for  angles  between  90  and  180  degrees),  however,  the  MoM-PO  prediction  is  not  accurate  enough.  The 
discrepancy  between  the  pure  MoM  and  hybrid  MoM-PO  here  is  certainly  primarily  due  to  the  fact  that 
the  currents  over  the  parts  of  the  cylinder  surface  in  the  shadow  region,  which  are  quite  weak  as  compared 
to  the  currents  in  the  lit  region  but  contribute  significantly  to  the  (low-field)  radiation  in  the  back  region 
of  the  antenna  system,  are  set  to  be  zero  in  the  MoM-PO  simulation.  Using  the  two-fold  symmetry  of  the 
problem,  the  total  number  of  unknowns  for  the  antenna  system  amounts  to  6181  with  both  techniques.  In 
the  pure  MoM  analysis,  the  CPU  time  is  about  19  minutes  with  a  relatively  modest  PC  (AMD  XP-1700+ 
with  512  MB  of  RAM).  In  the  hybrid  MoM-PO  analysis,  the  number  of  unknowns  in  the  MoM  region  is 
20  and  that  in  the  PO  region  646 1 .  The  CPU  time  is  about  2  seconds  with  the  same  PC,  which  is  by  three 
orders  of  magnitude  faster  than  with  the  rigorous  (full  MoM)  higher  order  technique.  Note  also  that  the 
estimated  number  of  unknowns,  based  on  a  topological  analysis,  for  a  common  low-order  small-domain 
MoM-PO  solution  with  RWG  basis  functions  on  flat  triangular  patches  and  the  use  of  two-fold  symmetry 
is  more  than  65000  for  the  analysis  of  the  same  problem. 
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Fig.  5.  Higher-order,  large-domain  MoM-PO  Fig.  6.  Far  field  of  the  antenna  system  in  Fig.  5, 

geometrical  model  of  an  antenna  system  consisting  computed  by  the  full  MoM  and  hybrid  MoM-PO 

of  an  array  of  nine  A/2  dipoles  in  front  of  a  large  higher-order  techniques,  respectively,  in  the 

metallic  cylinder  of  finite  length.  horizontal  plane. 


5.  Conclusions 

This  paper  has  presented  higher-order,  large-domain  electromagnetic  simulations  based  on  the  finite 
element  method,  method  of  moments,  and  physical  optics.  The  implemented  analysis  tools  are  double¬ 
higher-order  computational  techniques  that  combine  higher  order  geometrical  modeling  and  higher  order 
field/current  modeling,  which  have  been  shown  to  be  equally  important  for  accurate  and  efficient  FEM, 
MoM,  and  PO  computations. 
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